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ABSTRACT 


Tht‘ current methodology for the aseismic design of structures is largely 
based on the concept of ensuring structural safety during a single earthquake 
event This event is usually the most severe event, during which availability of 
certain minimum ductility is ensured in the structure This design methodology 
does not consider the occurrence of several not-so-severe earthquakes and damage 
accumulation during each of these earthquakes, as may be necessary in the areas 
of moderate to lugh seismicity, and thus the actual life of the structure may 
sometimes be significantly less than the design life. In this study, a new approach 
heis been proposed for the estimation of the actual design life of a single-degree- 
of-freedom (SDOF) structure in a given seismic environment such that at the 
end of this life, the structure is expected to attain a specified level of damage 
This approach is based on the estimation of expected number of earthquakes 
at each som'ce with the help of Gutenberg-Richter relation and time-dependent 
hazard model, and on the estimation of structural damage during each of these 
events. The estimation of structural damage is based on the use of scaling models 
for Fourier spectrum, peak ground acceleration, and strong motion duration, on 
finding the properties of the equivalent linear SDOF oscillator, and on using 
the damage model proposed by Park and Ang, The structure is assumed to be 
hysteretic, elasto-perfectly plastic oscillator, and the strength and stiffness of this 
oscillator are assumed to degrade during each damaging event. A case study has 
been carried out for illustrating the proposed model, and it has been shown how 
this may be used to determine the design force levels for maximum allowable 
damage at the end of the design life. 
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CHAPTER I 
INTRODUCTION 


1.1 General Introduction 

The earthquake-resistant design of structures according to the existing 
design philosophy aims to ensure that during their lifetime, the structures re- 
sist the maximiun possible earthquake without collapse. The magnitude of this 
earthquake is estimated based on the available statistical data on past earth- 
quakes If the past data is not available, fault rupture parameters are used to 
assess this magnitude, or in other cases, it is based on the expert Judgement. 
The cpicontrnl distance of tlie fault at whk li tills earthqiiako occur.s is normally 
taken as tlie closest distance of this fault from the site. By using the estimated 
values of magnitude and epiceiitral distance with an attenuation relationship, 
peak ground acceleration (PGA) is estimated and the structure is designed for 
the level of forces consistent with this PGA This design approach based on the 
single eaithquake of largest magnitude is termed as the ‘scenario earthquake* ap- 
proach. This approach may be unsafe in those areas where several earthquakes 
of smaller magnitudes are also likely to occui' besides the event for wliich the 
structure has been de.signed Tlus may also bo unacceptable wlien the stiiicture 
sliould not suffer damage beyond a critical limit due to the heavy financial losses 
associated with interruptions in business activity. Due to the first possibility, 
the cumulative damage in the structure may become so large during the .smaller 
magnitude oartliquakos only that the structiue may become unusable before its 
design life is completed. 
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For reliable evaluation of the design levels wliicU will be consistent with 
the allowable damage at the end of tlie design life, it is necessnty lio aecoimt. 
for the seismicity of the area iu a coinpreheiisive manner and to estimate the 
nunibor of difh'icnt iiingnitudc oventa with proper apatial dlatribution aroiiiid 
the site. It is also necessary to estimate damage due to each of these events and 
to get cumulative damage at the end of design life of the structure for all of these 
events 


Seismicity represents the expected rate of occurrence of earthquakes with 
different magnitudes. Along a fault, this is generally obtained by studying tiic 
frequency-magnitude relationship using the available information on the past 
earthquakes (Gutenberg and Richter (1942)), or by fault rupture parameters. 

Early analytical methods to determine the seismic hazard at a site were 
based on PGA and the return periods were calculated with respect to this pa- 
rameter only. Housner (1969) studied the attenuation of PGA in several regions 
of the United States and pro.sented his results graphically in torms of fault length. 
Esteva and Villaverde (1973) derived expressions for PGA and peak ground ve- 
locity (PGV) on the basis of accelerations repotted by Hudson (1971, 1972a,b). 
Esteva (1974) gave a procedure to combine the geological, geophysical and all 
other non-statistical evidences for producing a mathematical model of seismicity 
of the given area Attempts have also been made by Crouse (1973), Hudson and 
Udwadia (1973) and Salt (1974) to predict analytically the characteristics of mo- 
tions on different soils. Yegulalp and Kuo (1974) applied the theory of extreme 
value for estimating the probabilities that given magnitude.? will be exceeded in 
given time intervals. 
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Classical methods of time-series analysis have been applied by many re- 
searchers (e g , see Knopoff (19G4), Aki (1963), Vere-Jones (1970) and 3hlien and 
Toksoz (1970)) attempting to devise analytical models for random earthquake 
sequences Most commonly applied stochastic models of seismicity assume that 
events of earthquake occurrence constitute Poisson process and that the magni- 
tudes are independent and identically distributed. Some drawbacks of this model 
become evident in the light of statistical formulation and the analysis of the phys- 
ical process involved. The Poisson assumption implies that the distribution of 
the waiting time to the next event is not modified by the knowledge of the time 
elapsed since the last event, while the physical models of gradually acciunulated 
and suddenly released energy call for a more general renewal i^roce-ss such that 
the expected time to the next event decreases with time (Esteva (1974)). Statis- 
tical data shows that Poisson assumption may be acceptable when dealing with 
the large shocks through out the world (Ben-Maneham (I960)) 

Statistical analysis of waiting times between earthquake does not favoiu: 
the adoption of the Poisson model or of other forms of renewal processes Some 
of these forms assume that the waiting times are mutually independent with log- 
normal or Gamma distributions (Shlicn and Toksoz (1970)). Thereforo, onc-step 
memory models or other semi-Markov models given, by Patwardhan et al, (1980) 
and Kiremidjian and Suzuki (1987) with a time-dependent hazard rate such as 
the Weibull or the lognormal, have been used to evaluate the probabilities of oc- 
currence Other models which have been developed (e.g., see Vere-Jones (1970)) 
arc of the trigger type, i.e., the overall process of earthquake generation is con- 
sidered as fJio supeijiosil.ion of number of timo scries, each having a difl'cient 
origin, whore the origin times arc the events of Poisson process. 



4 


Quantification of damage in structures due to earthqualce loading is fa- 
cilitated by the evaluation of damage index. The damage indices are defined 
locally for each element and are then integrated globally for the entire structure. 
Generally, the damage indices are considered to be dimensionless parameters and 
are intended to range between zero for an undamaged structure and one for fully 
damaged or collapsed structure, with an intermediate value giving a measure of 
the partial damage. 

The two earliest and .simplest damage indicators were ductility and in- 
terstory diift Despite the limitations they have with respect to the effect of 
cyclic loading, tiicy arc most widely recognized ns the damage paramcl.crs be- 
cause they are easy to calculate and their physical meaning is simple. Zhu et 
al. (1988) studied the effect of peak ground acceleration to velocity ratio, a/u, on 
ductility of inelastic systems. They discussed the significance of a/u ratio as the 
parameter for the ground motion characterization from a seismological point of 
view based on their damaging capacity. Stephens and Yao (1987) studied the cu- 
mulative damage in the structure using displacement ductility In this case, the 
damage index was used to assess two test structures which showed a moderate 
correlation with observed damage 

As an improvement over the ductility ratio, Banon et al (1981) pioposed 
flexural damage ratio which accounted for the stiffness and strength degradations 
that occur under cyclic loading. Roufaiel and Meyer (1987) improved upon the 
model of Banon et al. (1981) as it was not consistent with the test data. They 
defined flexural damage ratio as the increa.se in flexibility at maximum deforma- 
tion divided by the inciease in flexibility at failure. This parameter showed a 



good correlation with the residual strength and stiffness of the specimens tested 
in flexure. 


Accumulation of damage under cyclic loading is usually considered either 
by a low-cycle fatigue formulation in which damage is taken as a function of 
the accumulated plastic deformation or a term related to the hysteretic energy 
absorbed during the loading Fajfar (1992) proposed equivalent ductility factors 
based on damage due to low cycle fatigue, and used those for the construction 
of inelastic spectra Chung et al (1987, 1989) used the fatigue formulation for 
finding strength loss due to each cycle of loading. They calculated the damage 
index by using Miner’s rule slightly modified by some weights wliich wore based 
on the average stiffnesses for specimens cycled to failuic. Waiig and Shah (1987) 
developed a model based on the concept of accumulated dainago Their model 
is capable of predicting the hysteretic behaviour of reinforced concictc beam- 
column joints, especially tlie strength and sbiffne^ss deterioration and the damage 
state of tile member 

Energy absorption was used by Gosam et al. (1977) as a measure of dam- 
age. He calculated damage as a cumulative energy ratio based on the assumption 
that if peak force has dropped below 75% of the yield value, the remaining ca- 
pacity becomes negligible. Park and Ang (1985) propo.sed a model based on the 
premise that the total damage due to caiLhqiiakc motion is hncnr combination 
of the damage caused by maximum deformation and the absoibed hysteretic en- 
ergy. Tliis model has been the most widely used model primarily because it has 
been calibrated for several damaged structures during the past eai'thquakcs (see 
Park et al. (1987)). 
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Anderson and Bertero (1991) showed that large amplitude atcoleration 
may not always cause appreciable damage by driving the structure into the non- 
linear range In fact, many repetitive non-linear excursions of relatively smaller 
amplitudes may be more damaging during an earthquake of larger duration 
Jeong and Iwan (1988) have studied the effect of eartiiquake duration on the 
damage of the structures by using Miner’s rule. They calibrated the model 
against observed failures of reinforced concrete columns. Basu and Gupta (1995) 
proposed a piobabilistic model to estimate the damage to the structure with 
a given ductility by using the order statistics of the lugher order peaks. This 
model relates the damage to the entire response process, and not just to the 
largest response, Dasu and Gupta (1990) proposed a stochastic tcclmiciiie for 
developing damage-based inelastic spectra for aseismic design of the structures 
which can be idealized as SDOF oscillators. For this purpose, the oscillators have 
been approximated by equivalent linear oscillators. 

Fajfar and Gasper§ic (1996) proposed a non-linear method called N2 
method for the aseismic design of reinfoiced concrete buildings. Tins method is 
mainly applicable to the structures predominantly vibrating in the first mode. 

This study proposes a new approach for the estimation of the design 
life of a SDOF structure by assuming time-dependent hazard model with log- 
normal distribution for the return period, and by considering the mo(Jel of Park 
and Ang (1985) for the estimation of damage during a single event. The ground 
motion for each earthquake event has been characterized by using the scaling 
models developed by Trifunac and co-workers (Trifunac and Brady (1975, 1976), 
Trifunac and Lee (1985)). 
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1.2 Organization 

This work is organized in tliree chapters following ttiis chapter. 

In Chapter II, the proposed approach has been formulated for the eval- 
uation of the total damage expected to take place in n SDOF structure over a 
given design life while the structure is situated in an area of known seismicity 
The information on seismicity has been used to find the number of earthquakes in 
each magnitude range which arc likely to occur at a fault dming tlie design life of 
the structure. For a particular earthquake event, the ground motion characteri.s- 
tics at the site of the structure have been estimated by using the available .scaling 
relation.ships, and the corro.sponding structural damage has boon estimated by 
considering oscillator to bo hysterotii, in nature. The ciiinulativi’ damage dur- 
ing the design life has been estimated by consideiing tlie most critical sequence 
of earthquake events and by considering ‘event to event’ stiffness and strength 
degradation for the oscillator. 

In Chapter III, the proposed approach has been illustrated by considering 
a hypothetical area having four faults with different activity rates. The growth of 
damage in the .structure with its age has been studied, and a parametric study lias 
been carried out to highlight the interdependence of various parameters wliich 
govern tlie design life of a structure, for a given maximum damage, in case of 
oscillators of different initial periods, 

A brief summary and conclusions of this study have been presented in. 


Chapter IV. 



8 


CHAPTER II 

FORMULATION FOR DAMAGE-BASED LIFE PREDICTION 


2.1 Brief Overview 

We are interested in obtaining the life of a structure at a site located in 
an area of known seismicity for a permissible cumulative damage in the structure. 
The structure is expected to be subjected to certain numbers of earthquakes of 
different magnitudes m tliis period which are consistent with the seismicity of 
the area. The number of earthquake events of different magnitudes which aie 
expected to occur in the area over a specified period from various nearby faults 
may be estimated by using available Gutenberg-Richter relationship and tune- 
dependent hazard rate. The ground motion at the site under consideration during 
a particular event may be characterized for a level of confidence in form, of the 
PSDF Even though scaling equations are available for directly estimating this 
PSDF (De (1997)), it is also possible to use the available scaling equations on 
Fourier amplitude spectrum (TVifunac and Lee (1985)), strong motion duration 
(Trifunac and Brady (1975)), and PGA (Trifunac and Brady (1976)) for this 
purpose. Assuming the ground motion process to be a stationary process, its 
PSDF may be estimated from the knowledge of Fourier amplitude spectrum and 
strong motion duration (Bendat andPiersol (1986)) and then, to account for non- 
stationarity, this PSDF may be scaled up or down uniformly at all frequencies 
so as to correspond to the same expected PGA as tlic PGA estimated by the 
regression relationsliip (Shrikliande and Gupta (1997)). Since the structures are 
designed to undergo significant inelastic deformations during the most severe 



and moderately severe earthtiuakc.s, it is necessary to account for the non-linear 
behaviour of structure in estimating the damage during each event For this 
purpose, stochastic linearization technique may be conveniently used to find the 
equivalent linear system pioperties in case of SDOF systems. Damage caused 
by each event may be then estimated by using a suitable damage model and 
information on ordered peak amplitudes (see, e.g,, Basu and Gupta (1995)). By 
assuming a sequence of various events (as expected during the design life of the 
structure) which causes maximum cumulative damage, the design life may be 
predicted for a specified critical damage after an iterative process. 

2.2 Seismicity and Expected Earthquake Events 

For the calculation of cumulative damage due to various seismic events, 
it is necessary to know the number of events of different magnitudes occuring in a 
given time. For this, we need to know the rates of occurrences of the earthquakes 
of different magnitudes for each contributing source. 

Let the magnitude, Mk be representative of all the earthquakes with 
magnitudes lying in a small interval, {M^ — AM/2, Mk + AM/2), centered at 
Mk. Let El denote the epicentral distance of the Ith source from the site under 
‘ consideration. The average rate of occurrence per year, n./jfc, of magnitude, Mfc, 
earthquakes at the Jth source may be obtained by using the following relationship 
given by Gutenberg and Richter (1942), 

logu/fc = aj - biMk . (2.1) 

Here, aj and are the constants estimated from the data of past earthquake 
records or from the known slip rate along the Ith. source If all the events cor- 



10 


responding to a magnitude class at a source are assumed to follow a Poissonian 
sequence of occurrence, 1/njjt would represent the moan return period of this oc- 
currence for Mk magnitude earthquakes at the Ith source However, it may not 
be realistic to assume a umform rate of occurrence for the earthquake events over 
a given period of time. The likelihood of occurrence increases with the increase 
in the quiscent period, i.e., the period of strain build-up following an earth- 
quake event. Hence, the return period should be distributed based on how much 
time has elapsed since the last event occurred on the source under consideration. 
Such a distribution follows from the assumed distribution of the return period 
by conditioning on the fact that no event has ocemred till the instant of time 
under consideration. Thus, if f{t) represents the probability den.sity function of 
the return period, and F{t) represents the corresponding cumulativo probability 
distribution, the hazard rate, 


Mt) = (2 2) 

would represent the probability of the event occurring in the time interval, (i, t-\- 
dt) on the condition that no event has occurred in the preceding t years, However, 
since the Poissonian distribution is a memoryless distribution, its use in Eq, (2.2) 
will not be able to model the quiscent period. We can alternatively consider 
the return period to be lognormally distributed such that for a given magnitude 
interval, the median of the (lognormally distributed) return period is equal to the 
expected value of the exponentially distributed return period as in the Poissonian 
model (see Todorovska (1994)). Thus, the hazard rate in case of kth magnitude 


class and Ith source may be written as 





0* 


Oki 1 


-<I) 

Int ->/* 
Qh 

i) 


(2.3) 
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where, Xik and Cifc respectively are the mean and standard deviation of the ran- 
dom variable, InTjjt, where, denotes the return period for the fcth magnitude 
event on the ith source. Further, ^(.) and #( ) are the standard normal den- 
sity and cumulative distribution functions. Due to the constraint as assumed by 
Todorovska (1994) and since the median of T/jt is equal to Aj*; in Eq. (2.3) 
becomes equal to — In(njjt). Further, as suggested by Todorovska (1994), Cik ~ 
0 2 may be considered to be a reasonable estimate. It may be noted that the 
distribution of the return period, Tj^, is based on the quiscent period with respect 
to the event of magnitude, , only It does not preclude the possibility of the 
occurrence of the eaithquakes belonging to the other magnitude classes before 
the occurrence of the Mk magnitude earthquake. 

If there are no events during the time of Tq years since the last event, 
the number of occurrences of the Mk magmtude events on the Ith source in Y 
years is given by 


nikiTo -1- Y\To) = / hik{T)dT . 

Jto 


(2.4) 


2.3 Estimation of PSDF 

For estimating the damage due to all the events as predicted by Eq (2 4), 
it is necessary to characterize the ground motion m terms of PSDF of the ground 
acceleration process for each event. This study considers the estimation of PSDF 
by using the known scaling relationships for Fourier spectrum, strong motion du- 
ration, and PGA in terms of the parameters hke earthquake magnitude, epicen- 



tral distance and geologic site conditions For the scaling of Fourier spectrum, 
the following scaling relationship as proposed by Triftinac and Lee (1985) has 
been considered, 


12 


logio FS{T) = M + AU(A, M, T) + -h b2{T)s 

fe3(r)A 


+ 


100 


+ h{T)+h{T)M^ (2.5) 


Here, M is the earthquake magnitude, A is the representative distance from 
source to station, and s (= 0 for alluvium, 1 for intermediate, and 2 for rock) 
represents the site condition, for the desired combination of event and site The 
representative distance, A, is expressed in terms of the epicentral distance, R, 
focal depth, H, fault size, i?, and coherence radius, 5o, as 




Further, in Eq, (2 5), 5’s represent the coefficients determined from a regression 
analysis at each period, T, and Att{A,M,T) represents the attenuation function 
given by 


AU(A^ M, T) = Aq{T) logic ^ R< Rq 

= A(r) logic Ac- (i?-i?o)/200 R>R(, (2.7) 


Here, Rq is the epicentral distance from the earthquake source at which the 
surface waves start dominating the ground motion, Ac is the representative dis- 
tance corresponding to jRoj and Aq{T) is the empirically determined attenuation 
function 


For the scaling of strong motion duration, following relationship given 
by Tnfunac and Brady (1975) is proposed to be used, 


T, = -4.88s + 2.33M -|- 0 U9R . 


( 2 . 8 ) 



It may be mentioned that besides relating the Fourier spectrum amplitudes with 
the PSDF amplitudes, the strong motion duration plays a key role in determining 
the total number of cycles and thus the structural damage during the earthquake 
excitation (Basu and Gupta (1995)). 

Assuming the earthquake ground motion to be a stationary process, the 
PSDF corresponding to the Mfc magnitude event occurring at the Zth source is 
calculated at frequency, tu, as 

G|*(u.) = (2 9) 

where, Zik{to) and Tik respectively arc the expected Fourier .spectrum and stiong 
motion duration for magnitude, and epicentral distance, i?j, as estimated 
from Eqs, (2 5) and (2,8), 

To include the effect of non-stationarity of ground motion, (jfAi(w) is 
proposed to be scaled up or down so as to correspond to the same expected PGA 
as that estimated by a suitable scaling relationship in case of magnitude, Mfc, 
and epicentral distance, Ri. For this purpose, it is proposed to use the scaling 
relationship given by Trifunac and Brady (1976) as follows 

log PGA = M + log Ao(i7) - log a(M) . (2.10) 

Here, Ao(H) represents the attenuation of acceleration with epicentral distance 
and a{M) represents the magnitude-dependent scaling constant as determined 
by a regression analysis. 

2.4 Equivalent Linear Oscillator and n.e.sponae Slailstics 


Once the PSDF is determined for an event of magnitude, Mjt, occur- 
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ring at the Ith source, it is convenient to estimate the structural response by 
replacing the given non-linear oscillator by an equivalent linear oscillator. Once 
the properties of the equivalent oscillator are determined by applying a suitable 
linearization technique, the PSDF of its response may be obtained simply by 
multiplying the ground PSDF with the squared modulus of the response transfer 
function. The response PSDF may now be used to estimate the relevant response 
parameters for calculating damage. 

Let us consider a hysteretic, SDOF structural system with initial fre- 
quency, u)n, post-yield frequency of vibration, aw„, and viscous damping ratio, C 
The equation of motion of this oscillator, when subjected to the base excitation, 
Up{t), may be written as 

X + 2(u)nX-i- = -Ug(t) ( 2 . 11 ) 

where, x and x denote the mass acceleration and velocity respectively. Further, 
f{x,z) denotes the acceleration term corresponding to non-linear hysteresis mod- 
elled as the combination of linear and elasto-plastic system. This is expressed as 
(Suzuki and Minai{1987)) 

f{x, z) = aujl -b (1 - a)ujlz (2 12) 

where, z is obtained from the following differential equation, 

i = ® [1 - - Xy) - U{-x)U[-z - Xy)] (2.13) 

In Eq (2 13), 17( ) denotes the Heaviside step function and Xy represents tlie yield 
displacement for the oscillator, By statistical linearization (see, e.g , Roberts and 
Spanos (1989)), the above non-linear oscillator can bo, replaced by an equivalent 



linear oscillator with, damping ratio, and natural frequency, We These param- 
eters can be evaluated by averaging the response over one cycle as proposed by 
Caughey (1960), and thus, we have 


15 


ulg{ay){l - a) 


Cwn . . /;^(1 - a)ujn 




io^ 




1 — erf 


\/2tr„ 


with 


Xr 


<Ty^ 




= A^(^-A-lsm2A)exp(-A^dA 


(2.14) 

(2.15) 

(2.16) 
(2 17) 


and erf( . ) representing the error function In Eqs (2 16) and (2.17), we have 


and 


^rms — 


A=cos- 1-- 




rdo) 




(2 18) 


(2.19) 


Uo ((ai|-tu2)2 + (C,a;,)2] 

It may be noted that an iterative procedure is usually necessary for obtaining 

the values of We and (^e Thus, with an initial guess of <Ty^ the values of cUe and 
Ce are obtained by using Eqs. (2.14) and (2.15), and then Cy is updated by using 
Eq. (2.16). After obtaining the linearized properties, Ug and Qy the response 
PSDF is calculated as 


TTl / ^ 


( 2 . 20 ) 


From the response PSDF, J5j)t(ai), now, the expected amplitude of the ith 
order response peak, i e., i?[£C(,)], may be estimated by using the order statistics 
approach as proposed by Gupta and Trifimac (1988) Thus, we have 

/ OO 

7?P(,)(77)d77 , 

■OO 


( 2 . 21 ) 



w 


where 

TV! 

''<■> ~ W - ‘VO - 1)1 ‘^*’'**' *' “ 

is the probability density function of the ith order peak. In Eq (2.22), 

= ;k 

and 

^OO 

P(v)= / J7(u)du (2.24) 

Jr) 

respectively are the probability density and distribution functions of the peaks 
in the displacement respon.se process, and 
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(2.25) 


is the total expected number of peaks in this process. In Eq. (2 23), e is the 
bandwidth parameter defined as 


e ~ 


^0-^4 ■^2 ^ 

A0A4 


(2.26) 


where An is, in general, the nth moment of the PSDF, £'jjt(cj), and is defined by 


/'OO 

Art = / w"^?{jfe(w)da>, n 

Jo 


0 , 1 , 2 ,... 


(2.27) 


2.5 Damage Model 

For the N ordered peaks with amplitudes, E[x(i)], 
estimated in the previous section, the damage is proposed to be calculated as 
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per the model of Park and Ang (1985) This model consists of simple linear 
combination of normalized deformation and hysteretic dissipated energy as 


Ait = — + 

®ii ^y^ii 


(2.28) 


where, x,n is the maximum displacement that the equivalent linear SDOF sys- 
tem would be subjected to during the base excitation, Xu (= where ji is the 
available ductility) is the ultimate displacement of the system under monotonic 
loading, represents the effect of cyclic loading on structural damage, EH rep- 
resents the total energy dissipation in the structure during the excitation, and 
Qy is the yield strength of the structme. Following the modifications introduced 
in this model by Kunnath ct al. (1992), following modified form of Eq (2.28) has 


been considered in this study 


= (2,29) 

In Eq (2.29), the first term is a simple, pseudo-static displacement measure. 
It takes 110 account of cumulative damage which is accounted for solely by the 
second term. The advantage of this model is its simplicity and the fact that it 
has been calibrated against a significant amount of observed seismic dEimage. 


The parameter, ccm) in. Eq. (2.29) refers to the expected largest peak 
amplitude, in the absolute response process, |a:(i)i. It may be easily 

estimated by following the same procedure as explained for i?[®(i)] in the previous 
section. However, in this case, the probability density function, p{i]) and N are 
to be respectively replaced by p{ri) {i] > 0) and N p{i]) and N may be expressed 
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and 



(2.31) 


For reinforced concrete elements, /3 can be defined as a function of the 
value of shear and axial forces on the section and the total amount of longitudinal 
and confinement reinforcement. By means of regression analyses using large sets 
of cyclic test data, Park and Ang (1985) have suggested relationsliips for (5. 
Many researchers (e.g., see Kunnatli et al (1990), Stone and Taylor (1993)) have 
also proposed regression relations for finding /3. This study considers /3 = 0.1 
which has been considered as a default value in the software, IDARC (Kunnath 
et al. (1990)) 


Assuming the response to be a narrow-banded process, the total number 
of cycles in the response may be considered to be equal to the total number of 


positive zero crossings given by 


AT = 

271 Ao 


(2 32) 


However, since the actual response process is not exactly a narrow-banded pro- 
cess, the total number of estimated peak amplitudes will be more than Hq Hence, 
it is assumed that the largest TVq only of those will contribute to the energy dis- 
sipation, EH. Thus, for o = 0, the total oneigy tlis.sipation is calculated by the 


following equation, 


EH = 5^4Qy(i?['i:(i)] - Xy) . (2.33) 

t=l 

The structural damage during the event of magnitude, Mk\ occurring at the Rh 
somce now follows from Eq, (2.29) as 


aJv(/i - 1) 


(2.34) 
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It may be noted that the parameter, Xy or Qy, is to be assigned suitable 
values depending upon the available ductility in the system. It is the common 
practice in the aseisrnic design to provide lower yield strength, Qy, than what 
would be required for the elastic response during the most critical earthquake in 
the lifetime of the structure. Qy is obtained by an iterative procedure in such a 
way that the maximum displacement of the (non-linear) structural system under 
the critical excitation is yield displacement, Xy, times the available ductility It 
may be noted that due to greater damage a.ssociated with greater reduction in Qy 
(compared to the elastic strength), liigher value of Qy than that calculated for a 
given ductility may sometimes have to be chosen to limit damage to a specified 
level during the most critical event. 

2.6 Cumulative Damage and Life Prediction 

Once the expected damage for each event is estimated as explained in 
the previous section, the expected cumulative damage in the design life of the 
structure may be estimated by considering the occurrence of events of various 
magnitudes at each source as discussed in Section 2.2, end by adding the in- 
dividual damage due to each of these events. It may noted that there will be 
no damage due to those events for which the expected maximum displacements 
do not exceed the yield displacement, Xy. Since the proposed seismic hazard 
model does not give any information about the sequence of the expected events 
during the design life of the structure and sequence is important for a stiffness 
and strength degrading structure, following approximate procedure is suggested 
to calculate the progressive damage during the de.sigu life 
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The expected number of events for magnitude, at /tii soiuce are 
estimated for first 5, 10, 15, .. years FVom this, the expected number of events 
in a particular block of 5 years is estimated by subtracting the number till the 
previous block from the number estimated till this block. For example, for the 
block starting at the 15 years of age of the structure, the expected number of 
events estimated for first 15 years is subtracted from that estimated for the first 
20 years. Witliin each block, all the anticipated events aie arranged in the de- 
creasing order of the damage they would cause to the just-built, i.e., undamaged 
structure. This is done to maximize the estimated damage during the block un- 
der consideration. The cumulative damage at the end of each block is calculated 
by adding the damage due to all the events suppo.sed to have taken place till that 
much age of the structure. 

The present study considers the fact that after the occurrence of each 
damaging event, there will be degradation hi the stiffness and strength of the 
structure. Following the logic that such degradation should be more in case of 
more damage and since no model is available yet to predict this ‘event to event’ 
degradation, following degradation model is proposed for this study. 

It is assumed that the yield displacement, increases by Aaiy after 
each damaging event such that 

•where, ki denotes the initial stiffness of the system and denotes the stiff- 
ness of the system after the occurrence of the ith event. It may be noted that 
Eq (2.35) is based on the assumption that the maximum value of Axy is Xy in 
case of zero stiffness. Now, assuming that the stiffness degradation is an increa.s- 
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ing function of the maximum displacement beyond the yield displacement during 
an event, the stiffness, kt+iy may be expressed as 




1 - 






(2.3(>) 


L J 

where, the parameter, 7, controls the late of ‘event to event’ degradation. This 
has been chosen arbitrarily as 0.1, considering the degradation in the moment 
capacity of the reinforced conciete members as modelled by Reinhorn et al. (1992) 
for a particular event. It may be noted that the strength degradation has been 
assumed to be governed by the stiffness degradation as in Eq (2 36) and by the 
increase, Aa3y, in the yield displacement as in Eq. (2.35). 


The above procedure gives the estimate of expected cumulative damage 
at the end of the design life of the structure. By an iterative procedure, the design 
life of the structure can be easily predicted for a specified cumulative damage. 


central LiBRARIf 

t T., KANPUW 
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CHAPTER III 

ILLUSTRATION OF THE PROPOSED APPROACH 


3.1 Brief Overview 

The approach proposed in Chapter II for the estimation of progiossive 
increase in the structural damage with time and prediction of the design life of 
the structure for a specified damage level has been illustrated in this chapter. 
One hypothetical seismic environment with complete knowledge of the seismic 
activity along each fault has been considered for this purpose. 

3.2 Site Data 

The structure is assumed to be located in a hypothetical aiea consisting 
of four faults. Two faults are located at a distance of 30 km each from the site, 
wliile the other two faults are located at 40 and 50 km each. The values of ot for 
these faults are taken to be 3.28, 4.09, 3.77 and 3,09 respectively while 5/ has been 
assumed to be uniformly equal to 0.86 for all the four' faults (Todorovska (1994)) 
The focal depths of the sources at all faults are assumed to be uniformly equal 
to 5 km Tlie entire area is assiuned to have alluvium geologic vsite condition, 

3.3 Results and Discussion 

Figs, 3.1 to 3.3 respectively show the progressive increase in the expected 
damage, as indicated by the index, H, with the age of the structure for the 
structui'c design life, Y = 25, 50 and 75 years. It is assumed in obtaining these 



Age in years 


Figure 3,1 Progressive Damage for a SDOF Oscillator with 
1 sec Period and Y" = 25 years for ajj, = 2.0 and /a 
= 2,3 and 4. 





Damage Index, 
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Figure 3.2 Progressive Damage for a SDOF Oscillator with 
1 sec Period and Y = 50 years for aji, — 2.0 and p, 
= 2, 3 and 4. 


Damage Index, 



Figure 3.3 Progressive Damage for a SDOF Oscillator with 
1 sec Period and F = 75 years for = 2.0 and fj, 
= 2,3 and 4, 




2fi 


figuies that zero time has elapsed since the occurrence of last event on each fault, 
and that the total number of events expected dui'ing the design life are uniformly 
distributed over the vaiious lilocks of 5 yoais. Fmthoi, Llio stuicLuie has boon 
assumed to have 1.0 .sec; period, 5% damping, u = 0, and has been designed based 
on the largest magmtude earthciuake that Ls expected to occur during its do.sign 
life. The yield strength level has been fixed by taking q;/j — 2, where cxr denotes 
the ratio of the chosen yield strength level to the calculated yield strength level 
for wliich the maximum nou-liiiear displacement during the most critical event 
is design ductility, /a, times the yield displacement In each figure, three curves 
corresponding to /a = 2, 3 and 4 have been shown. For a given structure, the 
damage increases almost linearly with age which is consistent with the uniform 
temporal di.stributiou assumed for various events It may be ob.served further 
that with the increase in /a, the expected damage as indicated by the damage 
index, jD, increases in all the three cases Actually, with the increase m design 
ductility, the response reduction factors tend to become larger leading to greater 
reductions in obtaining the yield level of forces. In fact, in case of y = 25 years, 
no damage is expected for ^ = 2 at any age of the structure because here, the 
response reduction factor is approximately equal to aj^, and thus, the structure 
behaves elastically during all events due to high yield strength level Similarly, 
in Fig. 3.2, the structures with /a = 2 and 3 are not damaged at all for the first 5 
and 10 years of ago. Coming back to Fig. 3.1, it may observed that for /a = 3 and 
4, the structure reaches the critical damage level of 0.8 for collap,se (feee Ang et 
al. (1993)) well before the design life of 25 years. It reaches this level at the age 
of 18 and 12 years respectively. Thus, the structures with higher design ductility 
arc likely to be ciibioally damaged much before the end of design life. With the 
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increase ui the design UEe also, the possibility of premature collapse appears to 
increase as indicated by Figs 3.2 and 3 3 Even though the design levels are 
higher in case of the .structures with greater design life, the total damage is much 
higher for those due to the exposure to much greater number of damaging events. 
The,se trends will now be sliowii moie clearly through a paranietiic study. 

For the parametric .study, 28 SDOF oscillators witli the periods of O.OG, 
0 08, 0.10, 0.12, 0.14, 0.16, 0 18, 0.20, 0 25, 0 30, 0.35, 0.40, 0 45, 0 50, 0.60, 
0.70, 0.80, 0.90, 1.00, 1.20, 1 40, 1.60, 1.80, 2.00, 2 50, 3.00, 3 50, 4.00 have been 
considered, and the variations of cumulative damage, design life, de.sign ductility, 
and the fact, or, cv;j have Vioon .sl.udied with the variations in tlu* o.scillntor time 
period. In eacli case, two of the other parameters have been kept at their default 
values while the third parameter has been given different values. The default 
values of have been taken as: /i = 3, otn = 2, V = 50 years, and damage index 
= 0.8, while the different values given for obtaining diffeiciit curves arc: /i = 2, 
3 and 4; aji = 1.5, 2.0 and 2.5; Y = 25, 50 and 75 years; and damage index = 
0.4, 0.6 and 0.8. 

Figs. 3.4 to 3.6 show the variations in the design life respectively for 
vai’ying damage, design ductility, and ajt. Fig. 3.4 shows the possible design life 
values for differout values of ciitical damages. It is seen tliat the struct m'e should 
be designed for a lower design life if the damage is not allowed to cross a lower 
level It is also seen that for the allowable damage level of 0 8, the stiff o.scillators 
can be designed for much greater life than the flexible oscillators wliile for the 
allowable levels of 0 4 and 0.6, both flexible and stiff oscillators have compaiable 
design life values. It is so because the response reduction factors are lower for 
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Time Period (sec) 


Figure 3.4 Design Life of a Set of Oscillators for an Allowable 
Damage, D = 0.4, 0.6 and 0.8, with ajj = 2.0 and 
/X = 3. 
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Time Period (sec) 


Figure 3*5 Design Life of a Set of Oscillators for /i = 2, 3 and 
4, with an = 2.0 and D = 0.8. 
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Figure 3.G Design Life of a Set of Oscillators for aji ~ i.5, 
2.0 and 2.5, with /x = 3 and D = 0.8. 
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more stiff oscillators, and thus, those are designed for higher design force levels. 
Fig 3 5 shows that with greater design ductility values, shorter design life should 
be considered for the designing of the structure This is because of the rea.soa that 
by increasing ductility, we reduce the yield level of the structure which leads to 
greater number of inelastic excursions and thus to more damage Fig 3.6 shows 
that by choosing higher design strength levels, the structures can be designed for 
significantly greater design life values 

Figs. 3.7 to 3.9 show the variations in the design ductility (with time 
period) respectively foi varying design life, a/t, and maximum damage. Fig 3.7 
shows that for a given oscillator, we can afford to have liigher values of design 
ductility for the structures with smaller design life with the allowable damage 
remaining unchanged. Further, stiffer oscillators on average can have higher 
design ductility as compared to the flexible oscillatois. It is seen in Figs. 3.8 
and 3 9 that a higher value of aji is compatible with greater design ductility for 
the same damage and that for a lower value of limiting damage, a lower design 
ductility should bo considered. 

Figs. 3.10 to 3.12 show how damage index depends on oscillator time 
period for varying design ductility, design life aud an respectively. Consistent 
with the trends as observed above, damage increases with the natural peiiod 
of the oscillator. It is liigher in case of structures with greater design ductility 
and same design life, and for structures with greater design life and same design 
ductility. Similarly, damage is seen to decrease in Fig. 3.12 for increasing linear 
design level. 

Figs. 3.13 to 3.15 show the variation of the mimmum value of an for 
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Figure 3,7 Design Ductility of a Set of Oscillators for Design 
Life, Y — 25, 50 and 75 years, with aji ~ 2,0 and 
D = 0.8. 
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Figure 3, 


9 Design Ductility of a Set of Oscillators for an Al- 
lowable Damage, D ~ 0.4, 0.6 and 0.8, with Y 
— 50 years and oLfi =2.0. 
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Figure 3.11 Damage Index of a Set of Oscillators for Design 
life, Y = 25, 50 and 75 years, with /i = 3 and an 
= 2 . 
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varying damage index, design ductility and design life respectively. All the tlu'oe 
curves show the need to provide liighcr linear design levels for more flexible 
structural hystoma to have same damage in the same design life. Other trentls 
are same as nhsnrved above, i.e,, there is a need to provide higher linear design 
l<>v<'l i.hiui (hut letiuiu'd on (he IhvhIs i>l icspunse u'tluc tion lu< l.ui', Ini i) Uniidug 
damage to lower levels, ii) greater design ductility, and iii) greater design life 
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Figure 3.13 aji for a Set of Os 
D = 0.4, 0.6 and 
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Figure 3.15 aji for a Set of Oscillators with Design Life, Y 
25 , 50 and 75 years for D = 0.8 and /i = 3. 
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CHAPTER IV 
CONCLUSIONS 


A new approach has been proposed in ttiis study for the estimation of 
the design life of a SDOF structure which is situated in a seismic enviionment 
The structme is expected to attain a specified level of allowable damage at the 
end of its design life, The proposed founulation is litisi'd on i) the time-dependent 
hazard model for the earthquake occurrences, ii) use of scaling equations for the 
estimation of ground motion diuing a given earthquake event, iii) linearizing the 
hystcretic SDOF oscillator by using the method of stochastic linearization, iv) 
use of modified Park and Ang’s model and order statistics approach for calcu- 
lation of the structural damage during each event, v) considering most critical 
sequencing of the earthquake events in a block of five years, and on vi) assuming 
a stiffness and strength degradation model for accounting for the degradation in 
the properties of the structure from one damaging event to another. 

An illustrative study based on a hypothetical example has qualitatively 
led to the following important conclusions. It has been found that the ‘single 
event-based’ conventional method of design may be inappropriate for ensuring 
safety in those areas where, besides the most critical oartlKiuake, several cartli- 
quakes of milder intensity may also occur during the design life of the structure 
For the usually adopted levels of force reduction from the linear levels, these 
earthqualces may generate sufficiently strong ground motions at tlic site of a 
structure so as to drive its response to be inelastic. Depending upon tlie damage 


4.3 


levels coiisiclcred acceptable by the owners in view of the functional icquireiiieiits 
of the structures, the reductions in the linear response levels should thus be much 
smaller and be not governed by the ductility of the system alone 
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